Tetraquark bound states in a Bethe-Salpeter approach 
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We determine the mass of tetraquark bound states from a coupled system of covariant Bethe- 
Salpeter equations. Similar in spirit to the quark-diquark model of the nucleon, we approximate 
the full four-body equation for the tetraquark by a coupled set of two-body equations with meson 
and diquark constituents. These are calculated from their quark and gluon substructure using a 
phenomenologically well-established quark-gluon interaction. For the lightest scalar tetraquark we 
find a mass of the order of 400 MeV and a wave function dominated by the pion-pion constituents. 
Both results are in agreement with a meson molecule picture for the /o(600). Our results furthermore 
suggest the presence of a potentially narrow all-charm tetraquark in the mass region 5 — 6 GeV. 

PACS numbers: 14.40.Bc, 14.40.Rt, 12.38.Lg 



Introduction 

The nature of the light scalar ++ meson nonet is still 
an issue under debate. Experimentally, the lowest-lying 
states in this channel are very broad which makes it dif- 
ficult to identify their structure and properties. While 
other multiplets are firmly established to be qq mesons, 
the lightest scalar mesons <x, k, ao(980) and /o(980) dis- 
play signatures that signal a potential strong non-gg com- 
ponent. Once abandoned from the particle data book, 
these states were reintroduced only in the last decade due 
to new experiments such as KLOE in the e~e + — > 7r°7r°7 
channel [1] or BES in the J — > um~ir + channel [2]; see 
also [3] for a thorough compilation. 

Recently, several approaches to data analyses utilizing 
conventional and modified Roy equations have deduced a 
pole mass and width of the order of m a w 450-M280 MeV 
for the (t//o(600) [1, 4], with error bars below the five- 
percent level. Of course, these analysis do not answer the 
question on the nature of these states. In this respect, 
a number of arguments surfaced over the years which 
question the qq nature of light scalar mesons. In the non- 
relativistic quark model, scalar quarkonia are p waves, 
yet the ++ nonet members are comparatively light and 
the /o(600) lies even below the 1" isoscalar (s— wave) 
lo meson. The mass ordering within the ++ nonet is 
puzzling as well: the lowest-lying state is the isoscalar 
/o(600) instead of the isotriplet do (980), and the mass 
degeneracy of a (980) and / (980) is hard to reconcile 
with their different flavor content. Finally, their decay 
channels disagree with a qq picture: both ao(980) and 
/o(980) couple to KK although only the latter contains 
strangeness, and the broadness of the dominant decay 
channel Jo (600) — > ttt: remains unexplained. 

The unexpected behavior of these light scalar mesons 
can be resolved in a tetraquark assignment which was 
introduced long ago by Jaffe [5]. Here, the oddities men- 
tioned above are naturally explained by the flavor struc- 
ture of the qqqq nonet. The mass spectrum deduced 
from the tetraquark nonet is inverse to that of the qq 
nonet, thus explaining the low mass of the isoscalar. A 
++ tetraquark carries zero quark orbital angular mo- 



mentum [6] , in agreement with the expectation that such 
s-wave states should be light. The putative decay chan- 
nels of the tetraquark nonet agree with the observed ones, 
i.e., the coupling of and /o to KK and r]ir is caused 
by their strange-quark content, while the broadness of 
the /o(600) is a consequence of its OZI-superallowed de- 
cay into two pions. These phenomenological arguments 
based on the group structure of the tetraquark flavor 
nonet are backed by effective theory studies and large- 
N c arguments (see e.g. [7-12] and references therein) as 
well as recent lattice calculations [13-15] which suggest 
a strong qqqq-comporient in the lowest-lying ++ states. 

In this letter we present the first results for tetraquarks 
in a covariant continuum approach based on the cor- 
responding four-body equation for two quarks and two 
antiquarks [16]. We construct a suitable representation 
of this system in terms of mesons and diquark degrees 
of freedom. While we choose this approximation for 
the sole reason of its numerical simplicity compared to 
the full four-body equation, its feasibility is well mo- 
tivated by recent results in the baryon sector. There, 
the analogous quark-diquark approximation to the nu- 
cleon's three-body Faddeev equation works well on the 
five-percent level [17-19]. The relevance of diquark com- 
ponents in the systematics of hadron physics has also 
been emphasized in Refs. [20]. It is furthermore impor- 
tant to note, that the meson and diquark degrees of free- 
dom that appear in our approach are not fundamental; 
they are obtained dynamically from quark-(anti-)quark 
Bethe-Salpeter equations (BSEs) and a model for the 
quark-gluon interaction which has been very successful 
on a phenomenological level [21]. We therefore deter- 
mine the properties of tetraquarks from the fundamental 
quark and gluon degrees of freedom of QCD. 

The Bethe-Salpeter equation for tetraquarks 
We start with the full four-quark Green function which 
satisfies the relation 



G — Gq -\- Gq T Gn 



(1) 



Here, Go is the product of four dressed quark propagators 
and T represents the connected and fully amputated four- 








FIG. 1: Tetraquark BSE in the meson-meson/antidiquark-diquark picture. The hatched amplitudes involve diquark quantities; 
the remaining ones are of mesonic nature. Single (double, dashed) lines are dressed quark (diquark, meson) propagators. 



quark scattering matrix. In our symbolic notation the 
multiplications in Eq. (1) represent four-dimensional in- 
tegrations over the appropriate number of momenta and 
all indices are left implicit. Tetraquark bound states with 
mass M appear as poles in the scattering matrix 



T 



P 2 + M 2 



(2) 



and thereby define the tetraquark's covariant bound- 
state amplitude ^, with ^ being its charge conjugate 
and P the total momentum. 

The Green function satisfies a Dyson equation which 
relates it to the four-body interaction kernel K via 



G = G + G KG T = K + KG T . 



(3) 



Substituting Eq. (2) into (3) and projecting onto the sin- 
gular part, one arrives at a homogeneous bound-state 
equation for the tetraquark amplitude: 



* = KGnHf. 



(4) 



This equation can be solved once the four-body kernel 
K is known. The kernel contains 2PI, 3PI and 4PI con- 
tributions [16]; in the following we adopt the successful 
strategy of Ref. [18] and neglect the latter two. 

To simplify the notation, we suppress the quark propa- 
gators by replacing Gq 1 G -> G, T G -> T and K G -)• 
K, so that Eqs. (1), (3) and (4) become 

G=1 + T=1 + KG, T = K{\+T), * = A"*. 

The remaining part of the kernel that involves only two- 
body correlations can be written as the sum of three 
terms: 



(5) 



where a and a 1 denote qq, qq or qq pairs, so that aal 
is one of the three combinations (12)(34), (13)(24) or 
(14) (23). Therefore, K aa ' describes the component of 
the four-body kernel where all interactions are switched 
off except those within the pairs a and a' . 

The absence of residual color forces between widely 
separated clusters, potentially generated by the per- 
muted interactions in Eq. (5), amounts to the condition 



Here, G aa ' = 1 + T aa > is the four-body Green function 
obtained from the kernel K aa >, whose scattering matrix 
satisfies 



Taa' — K a a' (1 -|- T a , 



(7) 



and G a = 1 + T a , G a i = 1 + T a i are the two-body Green 
functions for the individual pairs a and a', with 



T a = K a (1 + T a ) , T a , = K a , (1 + T a ,) . 



(8) 



The separability of the Green function G aa ' imposes the 
following structure on the kernel K aa ' [16]: 



K„ 



K a + K a > ~ K a K a ' 



(9) 



where K a and K a > are now elementary qq, qq or qq ker- 
nels. The relations (7-9) yield the scattering matrix 



T a + T, 



(10) 



from which Eq. (6) can be readily verified. 

In principle, the tetraquark bound-state equation (4) 
with the kernel of Eqs. (5) and (9) can be solved with the 
techniques used in Ref. [18] for the three-body equation. 
In practice, however, this is a very demanding task in 
terms of computation power. While we strive to attack 
this problem in the future, for now we resort to a further, 
simplifying approximation in the spirit of the nucleon's 
Faddccv equation and its reduction to a quark-diquark 
picture [17]. In analogy to the three-body case, we define 
Faddccv amplitudes ^ aa / via 



(ii) 



Upon projecting Eq. (7) onto the bound-state amplitude 
"J, one obtains Faddeev-Yakubovsky type equations [22] 
for the amplitudes ^ a a r - 



*aa' = Taa' (* 



2 ' / bb' cc' 



(12) 



where T aa ' is constructed from two-body scattering ma- 
trices according to Eq. (10). 

Except for the omission of genuine three- and four- 
body correlations, Eq. (12) is still exact. Its reduction to 
a two-body problem proceeds by assuming that the two- 
body T— matrices are dominated by meson and diquark 
pole contributions, 



G a a ! — G a G a ' 



(6) 



T a {qi, 92, Q) = -T a (qi,Q) D a (Q) f a (q 2 ,Q) , (13) 
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and that the internal spin-momentum structure of the 
Faddeev amplitudes factorizes: 

*aa>(p,q,q',P)=r a (q,Q)D a (Q) 

(14) 

x T Q')D a ,(Q')$ aa ,(p,P). 

Without loss of generality we can assign the labels 1, 2 
to the quarks and 3, 4 to the antiquarks. Then, for 
ao! = (12)(34), r a and r a < describe diquark and antidi- 
quark bound-state amplitudes and D a , D a > their respec- 
tive propagators, whereas in the case of aa' — (13) (24) 
or (14) (23) the involved objects are of mesonic nature. 
Quantities with a bar indicate charge-conjugated ampli- 
tudes. P is the total tetraquark momentum and p the 
relative momentum between its respective constituents. 
The separated internal momenta q, q' correspond to 
the relative momenta of the (anti-)diquarks and mesons. 
Their total momenta Q, Q' can take arbitrary values and 
thus an off-shell description for the meson and diquark 
amplitudes is necessary. 

Combining Eqs. (12-14) and furthermore neglecting 
the single-interaction contributions T a and T a i from 
Eq. (10) yields a coupled diquark- antidiquark/meson- 
meson BSE which is depicted in Fig. 1. It describes an 
effective interaction between two mesons, or between a di- 
quark and an antidiquark, via quark exchange. We take 
into account the mesons and diquarks with lowest mass, 
i.e., the pscudoscalar- meson and scalar-diquark channels. 
Constituents with other quantum numbers are certainly 
possible; however, since the corresponding calculations 
are very expensive in terms of CPU time we postpone 
their inclusion to subsequent work. Moreover, we only in- 
vestigate the lowest-lying tetraquark with quantum num- 
ber J p = + . The resulting diquark-antidiquark and 
meson-meson contributions to the tetraquark amplitude, 

$n(p,P) := $(12X34), , 1E , 

®m(P,P) ■ = $(13)(24) = -$(14)(23) > 

are flavor and color singlets and Lorentz scalars. 

It is noteworthy that our framework, Fig. 1, does not 
permit a pure diquark-antidiquark state in isolation; it 
can only occur in combination with meson-meson in- 
teractions. On the other hand, both equations can be 
merged to a single mesonic equation where diquarks ap- 
pear only internally. Thus, one may view the result- 
ing tetraquark bound state as a meson molecule with 
diquark-antidiquark admixture to its kernel. We expect 
this diquark admixture to be especially important for 
tctraquarks with masses larger than the sum of their me- 
son constituents. 

Mesons and diquarks from quark and glue 
In order to solve the tetraquark BSE of Fig. 1, we need 
to determine the (on-shcll) masses and amplitudes of the 
meson and diquark building blocks as well as a suit- 
able continuation off their mass shells. For diquarks, 
this problem has been dealt with already within the 
quark-diquark approach to the baryon three-body prob- 
lem [17]. The corresponding techniques are well devel- 
oped and their reliability can be judged from the good 



agreement of nucleon and A masses in the quark-diquark 
picture with results from the corresponding three-body 
problem [18, 28]. We therefore adopt this framework also 
for the diquark and meson amplitudes that appear in our 
setup. The technical details of these types of calculations 
have been described in many works, see e.g. [21, 23, 24] 
for reviews, thus we only give a short summary here. 

The starting point is the Dyson-Schwinger equation for 
the dressed quark propagator, 

S ap(P) = Z 2 (¥ + m o) a f3 + J Kaa'P'P S a <p(q) , (16) 

9 

with wave-function renormalization constant Z2 and bare 
quark mass tuq. The exact interaction kernel K, aa 'i3'0 
contains the dressed gluon propagator as well as one 
bare and one dressed quark-gluon vertex. The Greek 
subscripts refer to color, flavor and Dirac structure. In 
the rainbow-ladder approximation that we adopt here the 
kernel can be written as 

2 47ra(fc 2 ) rpjiv a v / 17 \ 

K-aoc'W = ^2 p I k 7aa' 7/3/3', C 17 ) 

where T£ v = — k^k v /k 2 is a transverse projector 
with respect to the gluon momentum k. Eq. (17) de- 
scribes an iterated dressed-gluon exchange between quark 
and antiquark that retains only the vector part ~ 7^ of 
the quark-gluon vertex. Its non-perturbative dressing, 
together with the one for the gluon propagator, is ab- 
sorbed into an effective coupling a(k 2 ) which is taken 
from Ref. [18, 25]. 

Chiral symmetry and the associated axial-vector 
Ward-Takahashi identity demand the kernel /C to appear 
in the corresponding BSEs as well. The meson BSE is 
given by 

T aP (p,P)= J K aoi ,p, {S{q + )Y{q,P)S{q-)} a , pl , (18) 
1 

with the Bcthc-Salpeter amplitude T(p,P) depending on 
total and relative momenta of the quark and anti-quark 
constituents, and q± = q±P/2. The corresponding equa- 
tion for diquarks is obtained by the substitution of an 
antiquark with a quark leg. Both meson and diquark 
amplitudes contain four different Dirac structures which 
we compute explicitly. 

With a given effective coupling, one determines the 
quark propagator in the complex momentum plane and 
subsequently solves the meson and diquark BSEs. Ac- 
cording to the techniques developed in Rcfs. [17, 24], the 
resulting onshell wave functions are analytically contin- 
ued to offshell momenta and the effective meson and di- 
quark propagators are computed from their T— matrix 
relations. Finally, all building blocks are put together in 
the tetraquark BSE. Our numerical techniques used to 
solve this BSE are only slightly non-standard and will be 
described in detail elsewhere. 
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FIG. 2: Mass of the tetraquark as a function of the pseudoscalar-meson mass (left) and the quark mass (right panel). 



Results and discussion 

Our result for the mass of the up/down ++ tetraquark 
state as a function of the pseudoscalar-meson mass is 
shown in the left panel of Fig. 2, together with a calcu- 
lation that only includes the meson-molecule component 
of the tetraquark. Clearly, the meson-meson component 
dominates. Except for the chiral region, the overall de- 
pendence of the tetraquark mass upon the pseudoscalar- 
meson mass is linear within numerical errors, as can be 
seen from the comparison with the linear fit included 
in the plot. The reason for this behavior is even more 
clear from the right panel of Fig. 2, where we plot the 
tetraquark mass as a function of the quark mass. The line 
represents a fit to the data including a constant, a square 
root and a linear term. Apart from the constant term 
this is the typical behavior of a Goldstone boson. The 
tetraquark thus grossly inherits the mass behavior of its 
dominating pion-molecule constituents, with deviations 
generated from their interactions via quark exchange. 

One of the main results of our present work is the value 
for the u/d tetraquark at the physical point, i.e. the left- 
most points in Fig. 2, where to ps = m v . We obtain: 



u/d 
X Tetra 



c (0++) = 403 MeV, 



(19) 



with an estimated numerical error of ten percent. This 
value is only somewhat lower than the real part of the 
mass of the tr// (600), m a w 450-M280 MeV determined 
recently from experiment using Roy equations [1, 4]. Our 
value for the mass of the scalar tetraquark should also 
be compared with the corresponding one for an ordi- 
nary quark-antiquark scalar bound state which may mix 
with the tetraquark components. In our rainbow-ladder 
approximation such a state has a mass of m q q(0 ++ ) = 
665 MeV. It is well known that corrections beyond 
rainbow-ladder increase this value into the 1 GeV range 
[26, 27], whereas the pion mass is protected. Since our 
tetraquark is dominated by its meson-molecule nature, 
we therefore expect it to be stable against corrections 
beyond rainbow-ladder, while at the same time the mass 



splitting between the tetraquark and the quark-antiquark 
scalar will increase. Consequently, our results suggest to 
identify the physical lowest-lying scalar state to be domi- 
nated by a strong tetraquark component, which is in turn 
dominated by pion molecule contributions. Due to the 
(Pseudo-) Goldstone nature of the pion constituents, our 
result provides a ready and natural explanation for the 
small mass and the large decay width of the tr//o(600). 

In the strange quark region at about mQ uar k = 80 MeV 
we also observe an all-strange tetraquark bound state 
at roughly ?™ T ctraquark(0 ++ ) = 1-2 GeV. Certainly this 
state will mix with its pure ss counterpart as well as the 
lowest lying scalar glucball making an identification with 
/o(1500) or / (1710) not possible without further studies. 

It is furthermore interesting to speculate about the 
existence of an all-charm tetraquark state. Because of 
its flavor-structure in our meson-diquark picture, such 
a state would be a mixture between a meson and an 
axialvcctor-diquark component. Since already the scalar- 
diquark contribution is very small, we expect the axi- 
alvector component to be completely suppressed due to 
its larger mass. This leaves only the dominant meson- 
molecule part. In Fig. 2 the largest pseudoscalar-meson 
mass corresponds to a quark mass in the charm region. 
We therefore read off the mass of an all-charm scalar 
tetraquark state to be at 



Tetraquark 



(0 



5.3 ± (0.5) GeV , 



(20) 



where the error is a guess based on our numerical and 
systematic uncertainties. This mass is considerably lower 
than the 6.2 GeV obtained in simple model calculations 
[29, 30]. It is also much lower than the r/ c threshold. Po- 
tential decay channels into D mesons and pairs of light 
mesons necessarily involve internal gluon lines. The re- 
sulting decay width may therefore be rather small. 

Further results for tetraquark states with unequal mass 
constituents are numerically more demanding than the 
ones presented here and will only be available for a future 
publication. 
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